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develop  an  understanding  of  the  bit-error  rate  to  SNR  relationship  for  a  multiple-user 
CSK  optical  signal.  Using  the  theoretical  performance  as  a  guide,  a  MATLAB  model  was 
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I.  INTRODUCTION 


The  majority  of  data  transported  through  modem  U.S.  Naval  platforms  is  carried 
on  copper  wire  or  coaxial  cable.  As  the  amount  of  data  being  transferred  dramatically  in¬ 
creases,  the  need  for  a  more  efficient  methods  of  transport  arises.  Given  the  tremendous 
increase  in  bandwidth  of  fiber  optics,  its  use  to  transport  data  in  naval  platforms  is  be¬ 
coming  standard.  Not  only  does  fiber  optic  communication  offer  a  tremendous  band¬ 
width,  it  offers  many  other  desirable  characteristics  that  include  resistance  to  electro¬ 
magnetic  interference  (EMI)  and  an  excellent  signal-to-noise  ratio  that  make  its  use  on 
naval  platforms  extremely  attractive. 

To  design  the  information  flow  in  naval  vessels  of  the  future,  it  is  important  to 
study  methods  of  optical  data  transmission  and  explore  ways  to  improve  upon  and  exploit 
the  tremendous  bandwidth  available.  In  this  thesis  an  alternative  to  Wave  Division  Multi¬ 
plexing  (WDM)  will  be  developed  and  analyzed.  This  alternative  relies  upon  the  use  of  a 
specialized  coupler,  the  Mach-Zender  coupler,  to  demodulate  the  multiplexed  optical  sig¬ 
nal  into  individual  signal  components. 

The  Mach-Zender  coupler  is  composed  of  two  2x2  optical  couplers,  a  delay  line 
and  a  phase  shift  element.  The  signal  enters  the  first  coupler  input  port  (port  1)  as  shown 
in  Figure  1.1  and  is  equally  divided  between  the  output  ports  A  and  B.  Signal  al  is  then 
delayed  one  bit  period  (T)  and  becomes  signal  a\ .  Both  signals  then  enter  the  second 
coupler  are  combined  and  then  equally  divided  to  the  output  ports  A  and  B.  The  signal 
b'2  is  then  phase  shifted  -  nil  radians  giving  the  outputs  of  the  Mach-Zender  coupler  of 
a2  and  b2 . 

Utilizing  the  outputs  of  the  Mach-Zender  coupler,  a  receiver  can  be  constructed  to 
demodulate  a  single  optical  Differential  Phase  Shift  Keyed  (DPSK)  signal  as  shown  in 
Figure  1.2.  The  outputs  of  the  Mach-Zender  coupler  are  followed  by  a  photodetector 
(PD)  to  convert  the  optical  signal  to  an  electrical  current  by  an  amount  proportional  to  the 
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Figure  1.1  The  Mach-Zender  Coupler 


responsivity  (91 )  of  the  photodiode.  This  current  signal  is  further  amplified  by  a  transim- 
pedence  amplifier  of  gain  Zq  and  input  to  an  integrator  over  one  period  (7).  The  energy 
of  the  two  channels  is  then  compared  and  a  symbol  (“1”  or  “0”)  is  chosen  based  on  the 
channel  with  the  greatest  energy. 


Figure  1.2  Demodulation  of  DPSK  Signal  Using  Mach-Zender  Coupler 

By  combining  multiple  couplers  with  varying  delays  (e.g.  T  and  27),  a  Code  Shift  Keyed 
(CSK)  signal  employing  Hadamard-Walsh  coding  could  be  demodulated.  This  orthogo¬ 
nally  coded  signal  would  allow  multiple  stations  to  use  a  single  fiber  optic  line  simulta¬ 
neously  as  is  currently  being  done  using  Wavelength  Division  Multiplexing  (WDM). 
With  WDM,  however,  there  is  a  significant  cost  in  creating  very  sensitive  filters  neces¬ 
sary  to  demodulate  the  signal  [Ref.  1].  Figure  1.3  shows  a  system  using  multiple  Mach- 
Zender  couplers  that  is  capable  of  decoding  four  Hadamard-Walsh  coded  signals. 
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By  increasing  the  number  of  Mach-Zender  coupler  stages,  this  receiver  can  be  con¬ 
structed  to  demultiplex  increasingly  complex  Hadamard- Walsh  coded  signals. 

In  Chapter  II,  the  theory  behind  the  Mach-Zender  coupler  will  be  discussed,  and 
from  this  theory,  a  mathematical  model  of  the  coupler  will  be  developed.  This  model  will 
include  both  time  and  frequency  domains.  Chapter  III  will  use  the  Mach-Zender  mathe¬ 
matical  model  of  Chapter  II  to  develop  a  single  user  DPSK  encoded  optical  receiver. 
From  this  DPSK  receiver  a  probability  of  bit-error  will  be  determined  based  on  the  sig- 
nal-to-noise  ratio  at  the  receiver.  Chapter  IV  will  then  propose  a  multi-user  Code  Shift 
Keyed  (CSK)  optical  receiver  and  develop  the  bit  error  rate  for  the  multi-user  receiver. 
Chapter  V  will  then  compare  the  theoretical  results  of  the  bit  error  rate  versus  signal-to- 
noise  ratio  in  the  single  user  receiver  with  those  obtained  by  a  computer  model  of  the  re¬ 
ceiver.  Chapter  V  will  continue  with  the  computer  model  to  conduct  a  sensitivity  analysis 
of  the  single  user  receiver.  Finally,  Chapter  VI  will  discuss  the  conclusion  of  this  thesis 
and  present  some  areas  for  future  research  with  these  optical  receivers. 
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II.  DEVELOPMENT  OF  A  MATHEMATICAL  MODEL  FOR  THE  MACH- 

ZENDER  COUPLER 

A.  OVERVIEW  OF  MACH-ZENDER  COUPLER 

In  the  following  analysis,  a  mathematical  description  of  the  Mach-Zender  coupler, 
as  has  been  previously  developed  at  the  Naval  Postgraduate  School  [Ref.  2],  will  be  re¬ 
viewed.  This  model  will  be  used  to  analyze  the  demodulation  of  an  optical  DPSK  signal 
using  a  single  Mach-Zender  coupler. 

As  previously  described,  the  Mach-Zender  Coupler  consists  of  two  2x2  optical 
couplers,  a  delay  line  of  some  delay  nT  where  n  is  an  integer  and  T  is  the  period  of  the 
signal,  and  a  device  to  provide  a  phase  shift  of  -  n  12  radians.  Port  1  of  the  first  coupler 
is  the  input  and  Port  A  of  the  second  coupler  and  the  output  port  of  the  final  phase  shift 
element  are  the  outputs. 

The  delay  of  the  Mach-Zender  coupler  can  be  accomplished  by  inserting  a  length 
of  fiber  equivalent  to  the  desired  delay.  Heinbaugh  [Ref.  2]  proposed  using  a  2x2  coupler 
with  a  single  input  and  a  single  output  to  obtain  a  n/2  phase  shift.  This  mathematical 
model  will  use  a  -n/2  phase  shift. 

The  input  signal  is  defined  as  s(t)  in  the  time  domain  and  S(f)  is  its  Fourier  trans¬ 
form  in  the  frequency  domain.  As  the  review  in  this  chapter  will  show,  the  output  in  the 
time  domain  is  s(t-T)-s(t)  for  port  A  of  the  second  coupler  and  s(t-T)+s(t)  for  the  output 
of  the  phase  shift  element.  This  is  as  illustrated  in  Figure  2.1. 


Figure  2.1  Signal  Throughput  of  the  Mach-Zender  Coupler 
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B.  DEVELOPING  A  MATHEMATICAL  MODEL 


In  this  section  each  element  of  the  Mach-Zender  coupler  will  be  analyzed  mathe¬ 
matically  and  an  expression  relating  to  that  particular  element’s  effect  on  s(t)  will  be  de¬ 
veloped. 

1.  Coupler  Analysis 

The  optical  coupler,  Figure  2.2,  has  four  ports  labeled  1,  2,  A  and  B.  Each  port 
can  be  used  as  an  input  or  an  output . 


a,-> 

1  A 

^  “3 

<-bj 

b3^ 

2x2 

cL — ^ 

2  B 

< —  &4 

< — 

b4-» 

Figure  2.2  Optical  2x2  Coupler 

Let  the  frequency  domain  inputs  of  the  2x2  coupler  of  Figure  2.2  be  described  by 
the  following  vectors  a  and  b  : 


V 

X 

a  = 

a2 

h  - 

b2 

a3 

u  — 

A. 

b4_ 

(2.1) 


The  input  vector  a  and  the  output  vector  b  are  related  by  the  scattering  matrix  S 
defined  as  follows  [Ref.  3]: 
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(2.2) 


0  0  17- 
0  0  j  1 
5=  1  J  0  0 
J  1  0  0 

where 


b=Sa 


"5" 

'o 

0 

1 

7 

2i 

^2 

l 

0 

0 

7 

1 

«2 

£ 

=  V2 

1 

7 

0 

0 

a3 

A. 

J 

1 

0 

0 

_a4 

(2.3) 


(2.4) 


Equation  2.2  can  be  simplified  if  the  matrices  of  Equation  2.4  are  partitioned  as 
shown  below: 


'o 

0 

:  1 

7 

a, 

h2 

l 

0 

0 

:  7 

1 

a2 

bj 

=  V2 

T 

y 

7  o' 

'0 

y 

A. 

.j 

i 

:  0 

0 

A»_ 

(2.5) 


After  partitioning,  Equation  2.5  can  be  rewritten  as  follows: 


"!1_ 

1 

"1 

7'" 

A" 

=  V2 

.7 

1 

A. 

'63' 

1 

'1 

7 

A" 

A. 

_7 

1 

A. 

(a3  +  jaA )  /  \{2 
( ja3  +  a4 )  /  yfl 

(a,  +ja2)/yf2 
O'a,  +a2)  / 


(2.6a) 


(2.6b) 
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Observe  that  in  the  ideal  2x2  optical  coupler,  no  part  of  the  input  signal  is  re¬ 
flected  back  to  the  input  ports.  Specifically,  when  the  only  inputs  are  a,  and  a2 ,  the  out¬ 
puts  are  confined  to  b3  and  b4  as  in  Equation  2.6b.  That  is  b3  and  b4  consist  of  only  in- 
phase  and  n  /  2  phase-shifted  components  of  5,  and  a2 . 


Figure  2.3  Single  Input  2x2  Coupler 

To  apply  this  result  to  a  practical  use  of  the  Mach-Zender  coupler,  consider  Figure 
2.3  in  which  one  input  signal  is  applied  to  one  input  port,  and  the  remaining  port  is  termi¬ 
nated.  Equation  2.6b  now  reduces  to  become: 


'  ax  /  -s/2 

o']  /  a/2  1 

A. 

Jax  /  -s/2  _ 

1 

io 

J5* 

Ml 

i _ 

(2.7) 


We  see  that  the  input  a}  is  decreased  by  a  factor  of  l/ -s/2  at  ports  A  and  B  and  the  output 
at  port  B  is  shifted  in  phase  by  n  /  2  radians  from  the  input. 

2.  Delay  Element  Analysis 

As  shown  in  Figure  1.1  the  delay  element  represents  a  time  shift  of  T  in  the  time 
domain  which  corresponds  to  multiplication  by  eJ2vfr  in  the  frequency  domain.  The  input 
and  output  of  the  delay  element  can  be  modeled  in  the  time  domain  as  in  Figure  2.4  with 
the  output  a{  =  alej2>fr  in  the  frequency  domain. 
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a\ 


T 


Figure  2.4  Time  delay  element 


3.  Phase  Shift  Element  Analysis 

The  phase  shift  of  the  Mach-Zender  Coupler  can  be  accomplished  using  a  2x2 
coupler  with  only  one  input  and  one  output.  The  result  would  be  a  signal  shifted  by  n  /  2 
radians.  For  this  model,  a  phase  shift  of  -nil  radians  was  assumed.  The  phase  shift 
can  be  modeled  in  the  frequency  domain  as  in  Figure  2.5. 


a\ 


n 

2 


~  ja ,  -» 


Figure  2.5  Phase  Shift  Element 


C.  INPUT  SIGNAL  ANALYSIS 

Using  the  preceding  development,  the  function  of  the  Mach-Zender  coupler  can 
now  be  completely  described.  Consider  an  input  signal,  s(t),  applied  to  an  input  port  of 
the  Mach-Zender  coupler.  The  unused  input  port  (2)  is  terminated.  Let  S(f)  be  the  Fourier 
transform  of  s(t).  The  signal  can  be  applied  to  a  Mach-Zender  coupler  as  shown  in  Figure 
2.6. 


9 


Figure  2.6  Input  Signal  Through  Mach-Zender  Coupler 
Using  Equation  2.7  the  output  of  the  first  coupler  becomes 


V 

"S(/W 2~ 

A. 

_;s(/)/V2_ 

(2.8) 


After  the  delay  line  the  inputs  to  the  second  coupler  are  obtained  by  applying  the 
delay  illustrated  in  Figure  2.4.  The  inputs  to  the  second  coupler  are: 


~e-jW  o" 

"s(/)/V T 

~S(f)e-J2*IT  ItI2~ 

A'_ 

0  1_ 

jS(f)l47_ 

ys(/>/V 2 

(2.9) 


The  outputs  of  the  second  coupler,  b"  and  b\ ,  are  obtained  by  applying  Equation 
2.6b  with  b\  and  b4  as  the  inputs.  The  following  result  is  obtained: 


1 

'1 

f 

~ S(f)e-J2,!fT  /  4 2 

'l/2(S(/)e^-S(/))' 

A'. 

j 

1 

;s(/)/V 2 

JI2(S(f)el2*r -S<f» 

The  addition  of  the  -n  12  phase  shift  gives  the  final  result  in  the  frequency  do¬ 
main  as: 


"1  o ' 

~\n 

1 

~(S(f)e~J22fr  -  S(f)) 

A" 

_o  -j. 

_j !  2(S(f)e~J2nfr  -  $(/))_ 

~  2 

(SOX^  +  SCO). 
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By  taking  the  inverse  Fourier  transform  of  Equation  2.1 1  the  time  domain  result  is 
obtained: 


'my 

1 

s(t-T)~  s(t ) 

.m. 

~  2 

_s(t  -T)  +  s(t) 

(2.12) 


By  multiplying  the  matrices  used  to  represent  the  phase  shift,  delay  and  two  cou¬ 
plers  (Equation  2.13),  we  can  derive  a  general  equation  (Equation  2.14(a),  2.14(b))  where 
the  output  of  the  Mach-Zender  coupler  can  be  readily  determined  from  the  inputs  al  and 
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_K_ 

0  -j. 

V2 

J  1. 

0  l. 

V2 

J  lJL«2. 

Phase  2nd  delay  1st 

Shift  Coupler  Coupler 


X” 

_  i 

\e-j1^  - 1) 

(je-w +J)- 

"«i" 

X: 

“  2 

(e~j2’fr  + 1) 

_S2_ 

1 

2 


a,  ( e  J2,fr  - 1)  +  a2j(e  j2,fr  + 1) 
ax  ( e~J2,fr  + 1)  +  a2  (e'j2,fr  - 1) 


(2.14a) 


(2.14b) 


D.  SUMMARY 

In  this  chapter  the  optical  2x2  coupler  was  shown  to  combine  the  signals  of  its  in¬ 
put  ports,  split  the  signal,  and  send  the  result  to  the  output  port  such  that  both  signals 
leave  the  coupler  simultaneously.  It  was  determined  that  if  only  one  input  port  was  used 
then  no  part  of  the  input  signal  was  reflected  back.  It  was  also  shown  that  the  outputs 
from  the  Mach-Zender  coupler  are  (S(f)e~j2*fT  -S{f))/l  and  {S{f)e~J2nfr  +  S(f))  /  2 
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in  the  frequency  domain  or,  equivalently,  ( $(7  -  T)  -  s(t  j)  /  2  and  (s(t  -T)  +  s(t))  /  2  in 
the  time  domain. 

In  the  next  chapter  a  receiver  will  be  proposed  using  the  Mach-Zender  coupler 
that  is  capable  of  processing  a  DPSK  encoded  signal.  Using  the  results  of  Chapter  II,  the 
signal  will  be  traced  through  the  receiver  and  its  strength  relative  to  the  receiver  noise 
will  be  analyzed.  Finally,  a  probability  of  bit-error,  Pjj,  will  be  developed  based  on  these 
results. 
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III.  DPSK  DEMODULATION  ANALYSIS 

A.  INPUT  DPSK  SIGNAL 

In  Chapter  I,  an  optical  receiver  capable  of  decoding  a  DPSK  encoded  signal  was 
proposed  (Figure  3.1).  Chapter  II  then  developed  a  mathematical  model  (Equation  2.12) 
of  the  Mach-Zender  coupler  to  be  used  in  this  optical  receiver.  With  these  results,  the 
next  step  is  to  analyze  the  input  of  a  DPSK  signal  and  assess  the  receiver  performance 
versus  the  input  signal-to-noise  ratio  (SNR). 


Figure  3.1  Optical  DPSK  Receiver 

A  given  DPSK  signal  will  be  in  one  of  two  forms  (over  a  period  of  2T  where  T  is 
the  period  of  the  input  signal)  depending  on  the  information  bit  being  transmitted.  For  a 
bit  “1”  being  sent  where  0  <  /  <  2T ,  the  signal  will  be: 

Si(t)  =  ±ApT(t  -  T)cos(2xfct  +  6>0)  +  ApT(t)cos(2.7ifct  +  0O)  .  (3.1) 

For  a  bit  “0”  being  sent  where  0  <  t  <  2T ,  the  signal  will  be: 

So(t)  =  ±ApT(t  -  T)cos(2nfct  +  0O)  ±  ApT(t)cos(27fct  +  <90)  (3.2) 

where  fc  is  a  multiple  of  1/T  and  pT(t  -  iT)  with  i  an  integer  is  defined  as: 
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(3.3) 


(t  _  •t'i  _  1 1  W<t<(i+l)T 
Pt(  lT)  |0  otherwise. 

Due  to  the  symmetry  of  the  signal,  an  analysis  of  an  input  1  or  0  will  produce  the 
same  results.  To  simplify  the  calculations  an  input  of  0  (Equation  3.2)  will  be  considered. 


B.  RECEIVER  ANALYSIS 


1.  Output  of  Mach-Zender  Coupler 

In  Chapter  II,  the  output  of  the  Mach-Zender  Coupler  was  determined  in  Equation 
2.12.  Applying  Equation  2.12  to  the  input  signal,  given  the  transmission  of  a  0  (Equation 
3.2),  the  following  output  is  obtained  where  T  <  t  <  2T : 


~a2(t ) 

1 

0 

“  2 

_±  A  cos(2  7fct  +  0O)±A  cos  (2rfct  +  #0 ) 

Simplifying  this  expression  yields: 


~a2(t) 

0 

±  ^cos(2^,r  +  ^0) 

(3.5) 


which  are  the  outputs  of  the  Mach-Zender  coupler  given  an  optical  DPSK  input  signal. 


2.  Photodetector  Analysis 

The  next  stage  of  the  receiver  is  an  optical  photodetector  (PD)  which  converts  the 
optical  signal  into  an  electrical  current  signal.  The  PD  provides  an  output  current  equiva¬ 
lent  to  the  input  peak  power  (P=A^/2)  multiplied  by  the  responsivity  ( 91 )  plus  the  shot 
noise,  r|s(t)  [Ref.  4]. 

Taking  our  input  from  the  Mach-Zender  coupler  (Equation  3.5)  and  applying  it  to 
our  PD  gives  the  following  signal  in  the  receiver  over  the  interval  where  T  <  t  <  2T : 
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(3.6) 


"«3(0" 

%.( 0 

A(0_ 

n(A2/2)+rjSb(t)_ 

It  is  important  to  note  at  this  stage  that  the  shot  noise  produced  by  the  two  PD  de¬ 
vices  ( (/)  and  TjSi  (t) )  are  independent  of  one  another. 

3.  Amplifier  Analysis 

The  next  stage  of  the  receiver  consists  of  a  transimpedance  amplifier  with  gain  Zq. 
Along  with  a  gain  to  the  input  signal,  this  amplifier  will  introduce  additive  white  gaussian 
noise  (AWGN),  rjamp(t),  to  the  signal.  In  order  to  simplify  the  understanding  of  noise  in 
the  receiver  a  summation  of  total  noise  will  be  necessary. 

Both  the  amplifier  noise,  r\amp(t),  and  the  shot  noise,  77^  (0  and  rj  (t) ,  are  mod¬ 
eled  as  AWGN  with  zero  mean.  If  the  total  noise  in  each  individual  channel  is  rj(t)  then: 


=  Vamp  (0  +  r,s(t)  ■  (3.7) 

Therefore,  the  noise  content  of  each  channel  following  the  amplifier  stage  will  be  %  ( t ) 
and  T]b(t)  with  a  power  spectral  density  (PSD)  of  N().  Again,  the  noise  in  each  channel 
will  be  independent  of  the  noise  in  the  other  channel. 

Given  this  simplifying  assumption,  the  output  of  the  amplifier  is  simply  the  input 
(Equation  3.6)  plus  the  amplifier  noise  (  rjamp  ( t ) )  multiplied  by  the  amplifier  gain  (Zq): 


a4(0" 

ZovAO 

MO. 

_Z0 (9T42  /2  +  7]b(t ))_ 

(3.8) 


4.  Integrator  Analysis 

The  next  step  in  the  analysis  is  to  take  the  amplifier  output  (Equation  3.8) 
and  run  it  through  the  integrator.  This  will  give  us: 
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(3-9) 


a5(t) 
b5(t ) 


1 

\\ZMW 

0 

T 

J[(Z,9M2  / 2)  +  (Z0Tjb(t))]dt 


which  will  simplify  to: 


a5(t) 

b5(t) 


ZoyiA2T /2+  j[Z0Tjb(t)]dt 


0 


(3.10) 


To  continue  with  the  analysis  of  the  integrator  output  the  statistical  nature  of  the 
signal  at  this  point  will  be  examined.  Signal  a5  (t)  is  defined  as  follows: 

a5(t)  =  a5  +  Na(t)  (3.11) 

T 

where  a5  =0  and  Na  (t)  =  |[Z0^o(t)]<ir .  Since  Na(t)  consists  of  a  constant  multiplied 

o 

by  zero  mean  AWGN,  its  mean  is  also  zero.  Therefore  a5  represents  the  mean  and  Na  ( t ) 
represents  the  noise. 

Similarly  if  b5(t)  is  defined  as: 

b5(t)  =  b5  +  Nh(t)  (3.12) 

T 

where  b5  =  Z0ViA2T /  2  and  Nb(t)  =  J[Z07;A(r)]Jr ,  then  b5  represents  the  mean  and 

0 

Nh  (t)  represents  the  noise. 

Next,  for  a  gaussian  random  variable  such  as  our  zero  mean  AWGN  defined  by 
Na(t)  and  N}j(t),  the  variance  (  cr2 )  is  defined  by: 

°n,2  =  vNh2  =  Z02N0T .  (3.13) 
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Now  that  the  mean  and  variance  of  a5(t)  and  b5(t)  has  been  determined,  an 
analysis  of  the  decision  process  to  choose  the  signal  sent  can  be  undertaken. 

5.  Decision  Variable  Analysis 

With  the  mean  and  variance  of  the  signals  to  be  input  to  the  decision  element  de¬ 
termined,  an  examination  as  the  probability  of  an  error  in  signal  reception  can  begin.  A 
symbol  “0”  was  assumed  to  be  sent  to  our  device  initially,  therefore  an  error  will  have 
occurred  if  the  symbol  “1”  is  chosen  as  the  output.  This  probability  of  choosing  a  “1” 
when  a  “0”  is  sent,  Pr{l|0},  will  be  equivalent  to  the  magnitude  of  the  signal  in  channel 
a  being  greater  than  the  signal  in  channel  b  given  that  a  “0”  was  sent  or  Pr{ay>&5j0}. 

If  a  new  variable,  Y,  is  defined  where: 

Y  =  b5-a5  (3.14) 


then  it  is  apparent  that: 


Pr{a5  >  Z>5|o}  =  Pr{7  <  0|0} .  (3.15) 

From  the  analysis  of  the  integrator  output  it  was  determined  that  a5  and  b5  are 
both  independent  gaussian  random  variables,  and  therefore  Y,  as  the  difference  of  a5  and 
b5,  is  also  a  gaussian  random  variable.  Consequently  the  mean  of  F  or  F  can  be  ob¬ 
tained  as  follows: 


Y=b5-a5  =Z0ftA2T/2.  (3.16) 

Similarly  the  variance  of  Y ( aY2 )  can  also  be  obtained  as  follows: 

®r2  =o*.2  +  crNt2  =  Z02N0T  +  Z02N0T  =  2ZQ2N0T .  (3.17) 

Pr{l|0}  or  Pr{r  <  0|0}  can  be  defined  by  the  probability  density  function  of  Y  as 
follows: 
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0 

Pr{7<0|0}=  \fy(y\0)dy  (3.18) 

—oO 

where  the  probability  density  function,  fy  (y|0) ,  for  a  gaussian  random  variable  is  defined 
as: 


f,<y 10)  = 


no,, 


-exp 


•(( y-Y )2) 


lot 


(3.19) 


Substituting  Equation  3.19  into  Equation  3.18  gives: 


pr{r<o|o}=  )f,(yWy=  h=— «p 

-iV2^cr„ 


(( yzlll 

lot 


dy . 


(3.20) 


If  a  new  variable  x  is  defined  such  that  x  =  — —  then  dx  =  Using  this 


o„ 


new  variable  to  substitute  into  Equation  3.20  gives: 


Y !  (rv 


X 

T 


(  °’dx)~  4^) exp| 


r/c„ 


X 

T 


dx.  (3.21) 


Equation  3.21  is  now  in  the  same  format  as  the  Q  function,  Q(z) ,  where  z  =  — . 


o„ 


Therefore: 


Pr{l|0}  =  Q 


J  Z0WA2T/2 ^ 

=  Q 


Z'JiNJ 


(3.22) 
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C.  PROBABILITY  OF  BIT  ERROR 


The  probability  of  bit  error,  Ph ,  for  this  receiver  can  be  determined  as: 

Pb  =  Pr(l)  x  Pr{o|l}  +  Pr(0)  x  Pr{l|0} .  (3.23) 

If  a  “1”  or  a  “0”  is  equiprobable  (i.e.,Pr(l)  =  Pr(0) )  and,  if  by  symmetry  Pr{0|l} 

=  Pr{l|0} ,  then  substituting  Equation  3.22  into  Equation  3.23: 


(  1  2 

T  | 

2Nj 

(3.24) 


Further,  if  the  peak  power  (P)  is  defined  as  P=A2/2  then  Equation  3.24  will  simplify  to: 


=  Q\ 


5RP. 


Defining  the  voltage  signal-to-noise  ratio  (SNR)  as: 


SNR  =  —  =  ^9L42 
<x„  2 


(3.25) 


(3.26) 


then  Pb  becomes: 


Pb  =  Q(SNR)  .  (3.27) 

Figure  3.2  shows  this  dependence,  which  illustrates  the  rapid  decline  in  bit  error  rate  as 
the  signal-to-noise  ratio  exceeds  4  dB. 

D.  SUMMARY 

In  this  chapter  a  mathematical  model  for  a  single  user  DPSK  encoded  optical  re¬ 
ceiver  was  developed  using  the  analysis  of  the  Mach-Zender  coupler  from  Chapter  II. 
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Figure  3.2  Pb  vs.  SNR 

From  this  model  a  relationship  of  bit-error  rate  versus  signal-to-noise  ratio  was  then  de¬ 
veloped. 

In  Chapter  IV,  a  mathematical  model  of  a  multi-user  CSK  optical  receiver  will  be 
developed  in  a  very  similar  manner  to  the  development  in  this  chapter  for  a  DPSK  re¬ 
ceiver.  Again  as  in  this  chapter,  that  model  will  be  used  to  determine  the  performance 
(i.e.,  Pb  vs.  SNR)  of  the  CSK  multi-user  receiver. 
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IV.  CSK  DEMODULATION  ANALYSIS 


A.  INTRODUCTION 

The  preceding  chapter  developed  a  probability  of  bit  error  for  a  DPSK  demodu¬ 
lated  signal.  A  DPSK  demodulated  signal  implies  that  only  one  user  can  send  in  informa¬ 
tion  in  the  optical  fiber  at  a  time.  In  order  to  allow  multiple  users  to  utilize  the  fiber  si¬ 
multaneously,  the  signal  demodulation  must  be  modified.  As  was  introduced  in  Chapter  I, 
a  Code  Shift  Keyed  (CSK)  signal  employing  Hadamard-Walsh  coding  can  be  used  to  al¬ 
low  multiple  users  on  a  fiber  simultaneously. 

The  Hadamard-Walsh  code  is  a  straightforward  method  of  creating  multiple  or¬ 
thogonal  signals.  The  Hadamard-Walsh  basic  2x2  matrix  is: 


Observe  that  the  correlation  of  the  first  row  vector  with  the  second  row  vector  is  zero. 
The  matrix  of  Equation  4. 1  can  be  built  upon  to  create  additional  orthogonal  vectors  as 
follows: 


1  1 

1  -1 

1  1  = 

1  -1 


•1 


1  1 
1 

1  1 
1  -1 


(4.2) 


It  is  apparent  that  the  individual  elements  of  the  2x2  Hadamard-Walsh  matrix  are  used  to 
multiply  the  entire  2x2  matrix  and  Equation  4.2  can  be  written  as: 


1111 
1-1  1  -1 
11-1-1 
1-1-1  1 


(4.3) 
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Again  observe  that  each  row  vector  is  orthogonal  (i.e.,  has  a  zero  correlation)  to  all  oth¬ 
ers. 

The  4x4  matrix  of  Equation  4.3  consists  of  four  orthogonal  vectors  and  can  be  ap¬ 
plied  to  a  system  in  which  two  bits  or  two  users’  information  must  be  decided  simultane¬ 
ously.  It  is  this  matrix  that  will  be  used  to  design  an  optical  receiver  and  develop  the 
probability  of  bit  error  for  two  users  on  the  same  fiber  simultaneously. 

The  design  of  the  receiver  will  be  the  same  as  that  proposed  in  Chapter  I  and  is 
shown  in  Figure  4.1. 


B.  INPUT  CSK  HADAMARD-WALSH  ENCODED  SIGNAL 

A  given  two-user  CSK  Hadamard- Walsh  encoded  signal,  S(t),  will  be  comprised 
of  the  sum  of  two  orthogonal  signals,  one  each  from  channel  A  and  channel  B.  Therefore: 

S(t)  =  SA(t)  +  SB(t).  (4.4) 

S(t)  will  exist  over  a  bit  time  of  T=4T'  where  T'  is  the  period  of  the  input  signal.  Fur¬ 
ther: 


r  =  T-  (4.5) 

where  T  was  defined  in  chapter  three  as  the  period  of  a  DPSK  signal.  The  possible  signals 
from  channels  A  and  B  where  0  <  t  <  AT’  are: 

SlB(t)  =  +Apr  (t  -  37")  cos(2  7rfct  +  #0 )  +  Apr  ( t  -  2  7")  cos(2  jtfct  +  0Q ) 

+  Apr  ( t-T ')  cos  (2rfct  +  0O)  +  Apr  (t)  cos(2  jrfct  +  6'0)  (4.6) 

su(t)  =  +Apr  (t  -  3T')cos(2nfct  +  0O )  -  Apr  ( t  -  2  T’)  cos(2  tfct  +  0O) 

+  Apr  (i t-T ')  cos(l7ifct  +  0Q)~  Apr  ( t )  cos(2  jrfct  +  6>0 )  (4.7) 

S0B(t)  =  +Apr  (t  - 3T')cos(2nfct  +  0O)  +  Apr  (t  - 2 T') cos(2.7rfct  +  0b) 

-  Apr  (t-T')  cos(2  7rfct  +  0O)~  Apr  ( t )  cos(2  7rfct  +  0O )  (4.8) 
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Figure  4. 1  Multiple  Stage  Mach-Zender  Coupler  Receiver 


S0A  (t)  =  +Apr  (t  -  37”)  cos(2 Jtfct  +  0O )  -  Apr  it  -  2 7”)  cos(2 nfct  +  d0 ) 
-  Apr  ( t-T ')  cos(2nfct  +  00)  +  Apr  ( t )  cos(2 7rfct  + 0O ) 


(4.9) 


where  the  S  subscript  refers  to  the  bit  transmitted  and  the  channel  it  was  transmitted  from, 
respectively.  For  example.  Su\  (t)  refers  to  the  channel  A  component  of  S(t)  in  which  a  1 
was  transmitted  for  channel  A,  and  Sob  (t)  refers  to  the  channel  B  component  of  S(t)  for 
which  a  zero  was  transmitted  for  channel  B.  The  pr(t  —  iT')  term  is  defined  as  it  was  in 

Chapter  DI  (Equation  3.3).  The  frequency  fc  is  a  multiple  of  1/T'. 

Again  due  to  the  symmetry  of  the  signal  an  analysis  of  all  possible  inputs  is  not 
necessary,  as  each  will  produce  the  same  result.  Also,  as  each  signal  in  channel  A  is  or¬ 
thogonal  to  all  signals  in  channel  B,  there  will  be  no  interference  between  channel  A  and 
channel  B  signals  in  S(t).  To  aid  in  the  development  of  a  probability  of  bit  error,  Pb,  an 
input  of  a  1  in  channel  A  and  a  0  in  channel  B  will  be  considered.  This  input  corresponds 
to  Equations  4.6  and  4.7  respectively  and  the  input  signal  will  be: 

S(t)=S1A(t)+S0B(t).  (4.10) 


C.  RECEIVER  ANALYSIS 

1.  Output  of  First  Stage  Mach-Zender  Coupler 


In  this  multistage  receiver  the  first  stage  is  identical  to  that  analyzed  in  Chapter  III 
(Figure  3.1),  substituting  T'  for  T.  Again  as  in  Chapter  m  the  output  of  the  first  stage  is 
determined  by  applying  Equation  2.12  to  the  input  signal.  Equation  4.9.  Given  a  1  in 
channel  A  and  a  0  in  channel  B  the  following  first  stage  output  is  obtained  where 
T'<t<4T 

^(0  =J_COs(27 tft  +  60) 

|SB(0J  2 

A  detailed  derivation  of  this  output  is  contained  in  Appendix  A. 


■  2  Apr  ( t  -  37”)  +  2  Apr  (t  -  27j  -  2  Apr  (t-T') 
+  2Apr  (t  -  37”)  -  2Apr  (t  -  T') 


.  (4.11) 
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2.  Output  of  Second  Stage  Mach-Zender  Couplers 

The  second-stage  couplers  are  designed  with  a  2T'  delay  variable  and  conse¬ 
quently  the  analysis  of  their  outputs  is  different  than  that  of  the  first  stage.  However,  with 
only  a  slight  modification,  a  variation  of  Equation  2.12  can  be  used  to  describe  their  out¬ 
puts.  Specifically  for  a  delay  of  2T' Equation  2.12  can  be  rewritten  as  follows: 


a(t ) 

_  l 

~s(t  -  IT')  -  s(t) 

HO. 

~  2 

s(t-2T')  +  s(t)_ 

(4.12) 


Applying  Equation  4.12  to  the  input  of  the  second  stage,  Equation  4.11,  the  fol¬ 
lowing  output  is  obtained  where  3T'  <  t  <  4T' : 


's'o  At) 

0 

s;A(o 

_  i 

1  /  2(4  A  cos(2^fc  t  +  00)) 

So  AO 

~  2 

1  /  2(-4A  cos(2jtfct  +  60 )) 

_s;ao_ 

0 

which  can  be  simplified  to: 


's'oAO' 

0 

s;ao 

Acos(2jtfct +0O)) 

SoA  0 

-  Acos(2jtfct +90)) 

_s;ao_ 

0 

These  equations  represents  the  output  at  the  second  stage  Mach-Zender  couplers  given  an 
optical  multiplexed  CSK  Hadamard-Walsh  encoded  signal. 

2.  Photodetector  Analysis 

The  remainder  of  the  receiver  for  each  channel  is  identical  to  the  DPSK  receiver 
analyzed  in  Chapter  ID;  therefore,  the  analysis  will  be  similar.  In  the  next  stage  of  the  re¬ 
ceiver  the  optical  photodetector  (PD)  converts  the  optical  signal  into  an  electrical  signal. 
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The  PD  provides  an  output  current  equivalent  to  the  input  power  (P=A2/2)  multiplied  by 
the  responsivity,  ,  plus  the  shot  noise,  T|s(t). 

Taking  our  input  from  the  Mach-Zender  couplers  (Equation  4.13)  and  applying  it 
to  our  PD  gives  the  following  signals  in  the  receiver  over  the  interval  3 T'  <  t  <  AT' : 


X(o' 

V.JO 

Suit) 

9?(A!/2)+i),„(0 

Sob  it) 

S(A2/2)+I7s>,(() 

.Suit). 

ns„(» 

It  is  important  to  again  note  at  this  stage  that  the  shot  noise  produced  by  the  all  the 
PD  devices  are  independent  of  one  another. 

3.  Amplifier  Analysis 

The  next  stage  of  the  receiver,  as  was  the  case  in  Chapter  HI,  consists  of  a  tran¬ 
simpedance  amplifier  with  gain  Zq.  Along  with  a  gain  to  the  input  signal,  this  amplifier 
will  introduce  AWGN,  rfamp(t),  to  the  signal.  As  was  the  case  in  the  analysis  of  the  DPSK 
receiver  the  total  noise  in  each  channel  will  be  summed.  Since  both  the  amplifier  noise, 
rjamp(t),  and  the  shot  noise,  J]s(t),  are  AWGN  with  zero  mean,  the  total  noise  in  each  indi¬ 
vidual  channel  is  r)( t)  and  : 


ri(t)  =  'namp(t)  +  ris(t)  ■  (4-16) 

The  noises,  tJoaW,  liAft),  'HoB(t),  and  r\1B(t),  all  have  a  power  spectral  density  (PSD)  of  No 
and  the  noise  in  each  channel  is  independent  of  the  noise  in  the  other  channels. 

Given  this  simplifying  assumption  the  output  of  the  amplifier  is  simply  the  input 
(Equation  4.14)  plus  the  amplifier  noise,  T]amp(t),  multiplied  by  the  amplifier  gain,  Z0  and 
can  be  written  as: 
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(4.17) 


X'(0' 

ZqIJoa  i0 

S'uiO 

Z0(KA2/  2  +  ^(0) 

S'o'sit) 

Z0(MA2/2  +  ??0B(t)) 

.Slsit). 

Z0T}wit) 

4.  Integrator  Analysis 

The  next  step  in  the  analysis  is  to  take  the  amplifier  output  (Equation  4.17)  and 
run  it  through  the  integrator.  This  will  give  us: 


\[Z0Voa  (01  dt 
3  r 

47" 

\[z0m2  i2+rjXAmdt 

3  r 

AT 

j[Z0(9t^2  ll  +  VoBmdt 
3  r 

4r 

J[Z0»7i*(0]  dt 


which  will  simplify  to: 
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(4.18) 


(4.19) 


To  greatly  simplify  the  analysis  at  this  point,  notice  the  similarity  of  the  integrator 
output  defined  in  Equation  4.19  to  that  of  the  DPSK  receiver  integrator  output  defined  in 
Equation  3.10.  Specifically,  the  information  plus  noise  channels  in  both  cases  (S{"'(t) , 
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So'" (0 and  b5 (t) )  are  identical  as  are  the  noise-only  channels  ( S"B''(t),S""(t), and  a5(t)). 
Therefore,  the  statistical  nature  of  the  CSK  receiver  integrator  output  can  be  summarized 
as: 

Soa  =Sm+Nu  (4.20) 

Su=Sm+W«  (4-21) 

S„=Sm+Wm  (4-22> 

S,.=S,b+N,b  (4.23) 


where  the  mean  of  each  channel  can  be  described  as  S0A  and  Sw  (which  are  both  equal 
to  zero)  and  Su  =  S0B  =  Z0%42r'/2.  Further,  the  variances  of  each  channel  have  identi¬ 
cal  characteristics,  and  as  in  Chapter  HI,  each  channel  represents  a  gaussian  random  vari¬ 
able.  Therefore  the  variance  in  each  channel  can  be  defined  as: 

gn2=Z02N0T'.  (4.24) 

Now  that  the  mean  and  variance  for  each  channel  has  been  determined,  an  analy¬ 
sis  of  the  decision  process  can  be  undertaken. 


5.  Decision  Variable  Analysis 

With  the  determination  of  the  mean  and  variance  of  the  signals  to  be  input  to  the 
decision  element,  an  examination  of  the  probability  of  an  error  in  signal  reception  can  be¬ 
gin.  Due  to  the  symmetry  of  the  signal  the  probability  of  a  bit  error  in  channel  A,  Pb{A}, 
is  equal  to  the  probability  of  a  bit  error  in  channel  B,  Pb{B}.  Therefore  an  error  in  channel 
A  will  be  considered  and  the  result  will  be  the  same  for  channel  B.  If  a  symbol  1  was 
sent  in  channel  A,  an  error  will  occur  if  the  symbol  0  is  chosen.  This  probability  of 
choosing  a  0  when  a  1  is  sent  in  channel  A,  PrA  {o|l},  will  be  equivalent  to  the  magnitude 

of  the  signal  Sqa  being  larger  than  the  signal  Sja ■  A  new  variable  is  defined  as  follows: 
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(4.25) 


y=sia 


It  is  apparent  that  to  make  an  incorrect  decision  Y  must  be  less  than  zero  and  therefore: 

Pr,  >SU  |l}=  Pr,  (P  <0jl}.  (4.26) 


In  Chapter  HI  a  probability  of  bit-error  was  developed  (Equations  3.16  to  3.22) 
relating  one  information  plus  noise  channel  to  one  noise-only  channel.  Similarly,  in  this 
CSK  receiver,  channel  A  is  almost  identical  to  the  result  derived  in  Chapter  HI,  with  the 
only  difference  being  that  the  period  in  channel  A  is  now  T'  vice  T.  Using  the  result  ob¬ 
tained  in  Equation  3.24  and  substituting  T' for  T gives: 


Now  using  the  definition  of  T',  Equation  4.5,  and  rewriting  in  terms  of  T: 


pr,M=e 


iou2 


l4x2An 


(4.28) 


D.  PROBABILITY  OF  BIT  ERROR 

The  probability  of  bit  error  for  this  receiver  can  be  determined  as: 

P.=P,{4=P.{4  (4.29) 

Since  the  probability  of  an  error  in  channel  A  equals  the  probability  of  an  error  in  channel 
B,  Pb  can  be  generalized  by  considering  channel  A  only.  In  channel  A: 

Pfc{A}=  Pr,  (1)  x  Pr,  {0|l}+  Pr,  (0)  x  Pr,  {l|o}.  (4.30) 
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If  a  1  or  a  0  in  channel  A  is  equiprobable  (i.e.,  PrA(l)=PrA(0)=l/2)  and  if  by  symmetry 
PrA{0\l}=PrA{l\0}  then  by  substituting  Equation  4.28  into  Equation  4.30: 


J.9U*  —I— 
2  V4x2N0 


(4.31) 


If  the  power  (P)  is  defined  as  P=A2/2  then  Equation  4.31  will  simplify  to: 


Pb=Q 


1  T 


2Nn 


(4.32) 


and  if  the  signal-to-noise  ratio  (SNR)  is  defined  as  it  was  in  Equation  3.25,  then  Pb  be¬ 
comes: 


n=ef 


V 


(4.33) 


This  result  implies  a  3  dB  degradation  in  performance  when  compared  to  the  result  for 
DPSK,  Equation  3.26.  Figure  4.2  shows  a  plot  of  the  theoretical  performance  of.  the 
DPSK  and  the  CSK  receivers  and  at  bit-error  rates  below  about  3x1  O'3  the  DPSK  receiver 
is  degraded  by  3  dB. 


E.  SUMMARY 

This  chapter  built  upon  the  development  of  the  single-user  DPSK  receiver  from 
Chapter  IH  and  has  shown  that  the  performance  of  a  two-user  CSK  optical  receiver  is  3 
dB  less  than  that  of  a  single  user  DPSK  optical  receiver.  The  next  step  in  the  study  in  the 
use  of  the  Mach-Zender  coupler  in  an  optical  receiver  will  be  to  create  a  computer  model 
of  the  simple  DPSK  receiver  and  obtain  performance  data  from  the  model.  If  the  model 
performance  data  closely  resembles  the  theoretical  data  derived  in  Chapter  HI  then  vari¬ 
ous  components  of  the  model  can  be  altered  to  conduct  a  sensitivity  analysis  of  the  optical 
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SNR(dB) 


Figure  4.2  Pb  vs.  SNR  (theoretical  DPSK  and  CSK  receivers) 

DPSK  receiver.  Chapter  V  will  discuss  the  construction  of  the  computer  model  and  show 
the  results  of  the  sensitivity  analysis. 
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V.  SIMULATION  AND  SENSITIVITY  ANALYSIS 


A.  SIMULATION 


1.  Design  Approach 

The  third  chapter  developed  the  theoretical  probability  of  bit  error  for  a  DPSK 
optical  signal  demodulated  using  the  Mach-Zender  coupler.  If  a  simulation  of  the  optical 
and  electrical  receiver  components  can  be  constructed  using  MATLAB  code  and  the  sta¬ 
tistical  outcomes  of  the  simulated  receiver  at  various  signal-to-noise  ratios  accumulated, 
then  a  comparison  of  the  simulated  performance  to  the  theoretical  will  confirm  the  accu¬ 
racy  of  the  model. 


2.  Mach-Zender  Coupler  Simulation  Construction 

Previous  work  at  the  Naval  Postgraduate  School  [Ref.  2]  has  extensively  studied 
the  physical  characteristics  of  the  Mach-Zender  coupler.  Given  this  understanding  of  the 
coupler,  a  modularized  model  of  the  coupler  was  constructed  in  MATLAB  code.  Appen¬ 
dix  B  contains  a  copy  of  the  code  for  each  module.  Each  element  of  the  coupler  was  mod¬ 
eled  as  an  individual  function  with  varying  operational  characteristics. 

The  optical  2x2  couplers  were  designed  with  a  variable  input  scattering  matrix. 
This  allows  for  varying  amounts  of  coupling  of  the  input  signal(s)  into  the  output  signals. 
In  addition,  a  delay  element  was  constructed  capable  of  creating  any  desired  delay  onto  an 
input  signal. 

The  phase  shift  element  was  modeled  as  another  2x2  coupler  as  was  proposed  in 
the  phase  shift  element  analysis  of  Chapter  n.  A  drawback  of  using  another  2x2  coupler 
to  phase  shift  the  signal  is  the  3  dB  attenuation  of  the  signal  by  the  coupler.  In  order  to 
balance  the  final  output  power  of  both  channels  (£/  and  b"  of  Figure  2.6)  an  attenuation 

element  was  inserted  in  the  non-phase  shifted  output.  As  with  the  other  modules  the 
amount  of  attenuation  was  a  controlled  variable. 


3.  Electrical  Receiver  Element  Simulation 


The  next  stage  of  the  proposed  DPSK  optical  signal  receiver  converts  the  de¬ 
modulated  optical  output  of  the  Mach-Zender  coupler  into  an  electrical  signal  for  proc¬ 
essing  to  obtain  a  digital  output  signal  as  described  in  Chapter  m  Each  electrical  compo¬ 
nent  was  again  designed  as  a  module  with  varying  inputs. 

The  photodetector  can  apply  the  responsivity  (91 )  and  add  additive  white  gaussian 
noise  to  its  output.  The  transimpendance  amplifier  uniformly  amplifies  the  signal  and  also 
adds  additive  white  gaussian  noise.  The  integrator  integrates  the  signal  over  one  period. 
The  relative  strengths  of  the  signal  in  channel  A  and  channel  B  are  then  compared  to 
make  a  decision  of  a  “1”  or  a  “0”. 

4.  Statistical  Analysis 

The  MATLAB  model  of  the  receiver  was  repeatedly  run  at  various  SNR  to  obtain 
the  probability  of  bit  error  for  specific  SNR.  The  simulated  results  and  the  theoretical  re¬ 
sults  are  shown  in  Figure  5.1.  The  simulated  results  differed  by  about  4/10  of  a  dB  from 
the  theoretical.  This  indicates  that  the  MATLAB  model  was  accurate. 

With  the  accuracy  of  the  model  confirmed,  the  model  can  be  used  to  vary  the  op¬ 
eration  of  various  components  and  investigate  the  sensitivity  of  receiver  performance  to 
less  than  ideal  components.  Specifically,  the  amount  of  delay  and  the  phase  shift  were 
altered  to  non-ideal  values  and  new  receiver  probability  of  bit  errors  were  compiled. 

B.  SENSITIVITY  ANALYSIS 

1.  Delay  Error 

In  building  the  Mach-Zender  coupler  the  delay  element  was  implemented  by  in¬ 
serting  an  extra  length  of  fiber  equal  in  length  to  one  bit  period.  At  100  Mbps  Heinbaugh 
[Ref.  2]  showed  that  tolerances  of  ±0.5%  were  achievable.  Inserting  a  delay  error  of  of 
0.5%  into  the  MATLAB  model  showed  no  difference  in  Pb  vs.  SNR.  In  order  to 
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Pb  vs.  SNR  (ideal  receiver) 
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Figure  5.1  Pb  vs.  SNR  (ideal  receiver) 


acquire  measurable  data  delay,  errors  of  5%,  10%  and  20%  were  introduced  into  the 
model.  The  results  are  shown  in  Figure  5.2.  With  5%  delay  error  the  required  SNR  for  a 
specific  Pb  increased  one  half  of  a  dB.  A  10%  delay  error  resulted  in  over  a  3  dB  degra¬ 
dation  in  performance.  At  20%  delay  error  the  receiver  remained  a  random  number  gen¬ 
erator  (i.e,  Pb=0.5)  up  through  a  13  dB  SNR. 


Figure  5.2  Pb  vs.  SNR  (simulated  delay  error) 


1.  Phase  Shift  Error 

The  phase  shift  error  is  defined  by  the  ability  of  the  2x2  coupler  to  change  the 
phase  of  the  input  signal  while  keeping  attenuation  as  close  to  3  dB  as  possible.  In  the 
simulation,  attenuation  was  maintained  at  3  dB  but  the  effects  of  phase  delay  errors  of 
10°  and  30°  were  investigated.  The  results  are  shown  in  Figure  5.3. 


Figure  5.3  Pb  vs.  SNR  (phase  shift  error) 


For  a  10°  phase  shift  error  receiver  performance  degradation  was  minimal  and  almost 
zero.  For  a  30°  error  the  receiver  degradation  was  about  2  dB. 


C.  SUMMARY 

This  chapter  reviewed  the  construction  of  a  computer  model  in  MATLAB  code  to 
simulate  the  DPSK  optical  receiver.  The  performance  of  the  model  was  compared  to  that 
predicted  by  the  theoretical  results  obtained  in  Chapter  IH  and  the  model  performance 
closely  resembled  theoretical  performance.  Having  established  the  accuracy  of  the  model, 
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the  model  was  then  used  to  determine  the  sensitivity  of  the  DPSK  optical  receiver  given 
errors  in  delay  and  phase  shift  elements.  Chapter  VI  will  summarize  the  modeled  results 
and  discuss  areas  for  future  reasearch  regarding  the  use  of  the  Mach-Zender  coupler  in  an 
optical  receiver. 
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VI.  CONCLUSION 


A.  DISCUSSION  OF  RESULTS 

This  thesis  proposed  the  use  of  the  Mach-Zender  coupler  in  a  single-user  optical 
receiver  and  mathematically  illustrated  its  theoretical  implementation.  Next,  building  on 
this  implementation,  the  theory  supporting  the  use  of  the  Mach-Zender  coupler  to  de¬ 
modulate  multiplexed  fiber-optic  signals  was  developed  and  also  illustrated  mathemati¬ 
cally.  From  this  analysis,  we  created  a  computer  simulation  of  the  single-user  optical  re¬ 
ceiver.  Comparing  the  performance  of  the  simulation  with  the  expected  theoretical  per¬ 
formance  established  the  accuracy  of  the  model.  With  an  accurate  computer  simulation 
developed,  the  characteristics  of  the  model  were  then  altered  to  investigate  the  sensitivi¬ 
ties  of  the  optical  receiver  to  non-ideal  components. 

Chapter  IV  showed  theoretically  how  the  optical  receiver  could  be  expanded  to 
demodulate  a  two-user  multiplexed  fiber  optic  signal.  A  two-user  optical  receiver  experi¬ 
enced  a  3-dB  degradation  in  performance  as  compared  to  the  single-user  optical  receiver 
proposed  in  Chapter  in.  The  receiver  of  Chapter  IV  could  be  extended  to  increase  the 
number  of  users  multiplexed  on  a  single  optic-fiber  by  extending  the  number  of  Mach- 
Zender  coupler  stages  with  varying  delays. 

Chapter  V  discussed  the  receiver  simulation  using  a  MATLAB  model,  and  its  per¬ 
formance  closely  followed  theory.  This  indicated  that  the  model  was  accurate  and  that 
adjusting  parameters  of  the  model  to  more  closely  resemble  an  actual  Mach-Zender  cou¬ 
pler  should  give  reasonable  results  in  the  sensitivity  analysis. 

The  sensitivity  analysis  conducted  with  the  model  showed  only  one-half  dB  deg¬ 
radation  from  a  5%  delay  error.  As  was  discussed  in  the  Chapter  V,  a  5%  error  was  well 
above  the  0.5%  error  achieved  in  the  construction  of  the  Mach-Zender  coupler  by  Hein- 
baugh  [Ref.  2].  Therefore,  the  effects  of  delay  error  due  to  construction  on  receiver  per¬ 
formance  should  be  minimal. 

The  sensitivity  analysis  then  varied  the  phase  shift  error  up  to  30°  and  a  degrada¬ 
tion  of  less  than  2  dB  was  experienced.  A  phase  shift  error  of  30°  is  large;  larger  phase 
shift  errors  are  unlikely.  Therefore,  this  receiver  design  does  not  appear  to  be  sensitive  to 
probable  errors  in  phase  shift. 

This  sensitivity  analysis,  conducted  using  the  MATLAB  model,  indicates  that  this 
is  a  realizable  receiver  with  a  strong  potential  of  simultaneously  demodulating  multiple 
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signals  in  one  fiber  optic  line.  If  this  method  proves  to  be  feasible,  it  could  provide  a  cost- 
effective  alternative  to  demodulating  multiplexed  fiber-optic  signals  over  the  current 
technique  using  Wavelength  Division  Multiplexing. 

B.  FUTURE  RESEARCH 

Future  research  on  the  feasibility  of  using  the  Mach-Zender  coupler  should  focus 
on: 


1.  Using  the  existing  Mach-Zender  couplers  to  build  a  DPSK  optical-signal  re¬ 
ceiver  in  order  to  confirm  the  sensitivity  analysis  already  conducted  with  the 
MATLAB  model. 

2.  Constructing  a  CSK  optical-signal  receiver  to  investigate  the  reality  of  de¬ 
modulating  information  from  two  or  more  users  simultaneously  on  one  fiber. 
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APPENDIX  A:  ANALYSIS  OF  OUTPUT  OF  TWO  MACH-ZENDER 


COUPLERS  FOR  CSK  ENCODED  SIGNAL 


1.  SIGNAL  INPUTS 

The  following  two  signals  are  input  to  the  receiver  of  Figure  4.1  as  S(t): 

SlA  ( t )  =  +Apr  (t  -  3  T')  cos(2  jtfct +  60)~  Apr  (t  -  2  T')  cos(2 7tfct  +  60 ) 

+  Apr  (t  -  T')  cos(2 7tfct +0O)  -  Apr  ( t )  cos(2jtfct  +  60)  (A.  1) 

Sob  (0  =  +APr  ~  3T')  cos(2  jtfct  +  60)  +  Apr  (t  -  2 T')  cos(2 \jtfct  +  60 ) 

-  Apr  (t  -  T')  cos(2 jtfct +60)-  Apr  ( t )  cos(2  jtfct +60)  (A.2) 

due  to  the  orthogonality  of  the  signals  and  superposition  each  signal  can  be 
considered  separately  and  the  end  results  summed. 


2.  CHANNEL  A  INPUT  ANALYSIS 


First  consider  the  signal  for  channel  A.  To  aid  in  analysis,  cos(2 jtfct+60)  will  be 

factored  out  and  left  out  of  the  equation  until  the  end  result  is  obtained.  The  output  of 
the  first  coupler  with  a  delay  of  T'  for  the  channel  A  input  will  be  described  as  SA  (t) 
and  Sg  ( t ) .  The  superscript  refers  to  the  contributing  input  channel  and  the  subscript 
refers  to  the  output  channel.  Therefore  from  T'<t  <  AT'  the  output  is: 

SAA(t)  =  UslA(t-T')-S1A(t)}  = 

^[+Apr  (t  -  AT')  -  Apr  {t  -  3 T')  +  Apr  ( t  -  IT')  -  Apr  (t  -  T') 

-  [+Apr  (t  -  3 T')  -  APr  (t  -  IT)  +  Apr  (t  -  T')  -  Apr  (t)]] 

=  2  x  [—Apr  ( t  -  3 T')  +  Apr  (t  -  IT')  -  Apr  ( t  -  T')]  (A3) 
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SAB(t)  =  -[S1A(t-T')+S1A(t)\= 


1 


-  [+Apr  (t  -AT')-  Apr  (/  -  37')  +  Apr  (t  -  27')  -  Apr  (, t-T ') 


+  [+APr  ( t  -  37')  -  Apr  ( t  -  27')  +  Apr  (t  -  7')  -  Apr  (0]] 

=  ix2x[0]  =  0  (A.4) 


Equation  A.4  is  zero  and  can  be  disregarded.  Therefore,  applying  the  result  of 
Equation  A.3  to  the  input  of  the  second  stage  coupler  with  a  delay  of  2T'  gives  an  output 
from  37'  <t  <  AT'  of: 

SL'(t)  =  ^[SA(t-T)-SA(t)]  = 

j[-Apr  (t  -  57')  +  Apr  (t  -  47')  -  Apr  (t  -  37') 

-  \-Apr  (t  -  37')  +  Apr  (t  -  27')  -  Apr  (t  -  7')]] 
=  ~  x  2  x  [0]  =  0  (A.  5) 


S?A'(t)  =  ±[SA(t-T')+SA(t)]  = 

± [-Apr  ( t  -  57')  +  Apr  (i t  -  47')  -  Apr  (t  -  37') 

+  [-Apr  ( t  -  37')  +  ApT  (t  -  27')  -  Apr  ( t  -  7')] 

=  ^x2x  [-Apr  (t  -  37')]  (A.6) 

Equations  A.5  and  A.6  represent  the  output  of  the  second  stage  coupler  for  channel  A. 
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3.  CHANNEL  B  INPUT  ANALYSIS 


Next  consider  the  signal  for  channel  B.  To  aid  in  analysis,  cos(27tfct  +  60)  will  be 

factored  out  and  left  out  of  the  equation  until  the  end  result  is  obtained.  Therefore  at  the 
output  of  the  first  coupler  with  a  delay  of  7'  for  the  channel  B  input  from  7'  <t  <  AT'  is: 

SBA(t)  =  ^[SoB(t-T')-SoB(t)]  = 

~[+Apr  (t  -  AT')  +  Apr  (t  -  3 T')  -  Apr  (t  -  IT')  -  Apr  (t  -  T') 

2 

-  [+Apr  (t  -  3 T')  +  Apr  (t  -  IT')  -  Apr  ( t  -  T')  -  Apr  (?)]] 

=  ~  x  2  x  [-Apr  (t  -2T')]  (A.7) 


SBB(t)  =  ^[S1A(t-T')  +  SJA(t)]  = 

—[+Apr  (t  -  AT')  +  Apr  (t  -  3T')  -  Apr  (t  -  IT')  -  Apr  (t  -  T') 

+  [+Apr  (t  -  37')  +  Apr  {t  -  IT')  -  Apr  (t  -  T’)  -  Apr  (r)]] 

=  |  x  2  x  [+APr  (t  -  37')  -  Apr  ( t  -  7')]  (A.8) 

where  the  superscript  refers  to  the  contributing  input  channel  and  the  subscript  refers  to 
the  output  channel.  Equation  A.7  has  only  one  term  and  that  term  and  a  27'  delayed 
version  of  that  term  will  be  out  of  range  of  37'  <t<  AT'  and  can  be  disregarded. 
Equation  A.8  can  be  applied  to  the  input  of  the  second  stage  coupler  with  a  delay  of  2T’ 
gives  an  output  from  37'  <  t  <  AT'  of: 

SL\t)  =  ^[SB(t-T')-SB(t)}  = 

~[+ Apr  (t  -  57')  -  Apr  (t  -  37') 

jL 

-  [+Apr  (t  -  37')  —  Apr  (t  -  7')]] 

=  -^  x  2  x  [-Apr  (t  -  37')]  (A.9) 
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SlV(t)=\[SB(t-T')+SB(t)]  = 

~[+APr  (t  -  5T’)  -  Apr  (t  -  3T’) 

+  [+APr  (t  -  3T')  -  Apr  (t  -  r')]] 
=  ^  x  2  x  [0]  =  0 . 


(A- 10) 


Equations  A.9  and  A.  10  represent  the  output  of  the  second  stage  coupler  for  channel  B. 


4.  OUTPUTS 

The  output  of  the  first  coupler  in  equation  4.1 1  is  as  follows: 

'SA(tj 

MO. 

and  consequently  the  output  of  the  second  coupler  in  equation  4. 14  is  as  follows: 


Sj(t)  +  SB(0 
SA(t)  +  SB(t) 


X(  t) 

SoViO 

SuV) 

sfA\0 

s'oAt) 

SoViO 

MO. 

_C(0_ 

(A.  12) 


Substituting  Equations  A.  5,  A.6,  A.9,  A.  10  into  equation  A.  12  will  give  the  result  listed 
in  Equation  4.14. 
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APPENDIX  B:  MATLAB  SIMULATION  PROGRAMS 


This  appendix  contains  the  MATLAB  code  for  the  various  modules  of  the  DPSK 
optical  receiver  model. 


1.  OPTICAL  DPSK  RECEIVER 


%  set  input  paramters 


f  =1e6; 

T  =  1/f; 

Td  =  T; 

SM  =  [1,j;j,1]; 
L  =3; 

R  =1; 

No  =  1; 

Zo  =1; 

PWR  =  [8000]; 


%  frequency 
%  signal  period 
%  length  of  delay 
%  scattering  matrix 
%  attenuation  (dB) 
%  responsivity 
%  noise  PSD 
%  amplifier  gain 
%  power  levels 


t  =  linspace(0,T);  %  time  intervals  for  one  signal 


%  load  number  of  runs 
n=2000000; 


%  query  users  for  number  of  runs 

%n  =  input('enter  number  of  simulations  to  run»  '); 

for  p=1  :length(PWR) 

A  =  sqrt(PWR(p)); 

%  initialize  data  matrices 
xmtbit  =  zeros(1  ,n+1 ); 
rcvbit  =  zeros(1,n); 

%  initial  synchronization  signal  to  receiver 

S1=A*cos(2*pi*f.*t); 

phase  = 1 ; 

%  create  1  st  data  signal 
[S2,xmtbit(1)]  =  onechxmt(f, phase); 

S2  =  A*S2; 


%  update  phase  setting 
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if  xmtbit(l)  ==  1 
phase  = -1  "phase; 
end 

for  i  =  1:n 

%  create  1st  data  signal 
[S3,xmtbit(i+1)]  =  onechxmt(f,  phase); 
S3  =  A  *  S3; 

%  update  phase 
if  xmtbit(i+1)  ==  1 
phase  =  -1  "phase; 
end 

%  create  next  input  signal 
S  =[S3,S2,S1]; 

%  input  S  to  Mach-Zender  coupler 
[xa,xb]  =  mzcoup(S)T,Td,SM,L); 

%  input  MZ  output  to  RCVR 
[rcvbit(i)]  =  rcvr(xa,xb,R,Zo,No,T); 

%  progress  signal  to  next  period 

51  =  S2; 

52  =  S3; 


end 

%  resize  xmbit  to  remove  last  bit 
xmtbit  =  xmtbit  (1  :n); 

%  initialize  error  count 
errcnt  =  0; 

%  compute  total  nr  of  errors 
for  i  =  1:n 

errcnt  =  errcnt  +  xor(xmtbit(i),rcvbit(i)); 
end 

Pb(p)  =  (errcnt/n); 

SNR(p)  =  1 0*log1 0(R*PWR(p)*(sqrt(T/(2*No)))); 
end 
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%  procedure  to  save  data  to  .mat  file 
clear  P  S; 
load  ideal2dat; 

P=[P,Pb] 

S=[S,SNR] 

clear  L  A  R  No  Zo  SM  t  T  Td  i  n  Pb  SNR  errcnt  SI  S2  S3  rcvbit  xmtbit; 
clear  f  PWR  phase  p  xa  xb; 
save  ideal2dat; 
clear 

%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


2.  ONE  CHANNEL  TRANSMITTER 

function  [S, xmtbit]  =  onechxmt(f,phase) 

%  1  CHANNEL  DPSK  TRANSMITTER  one  bit  output 
%  [S, XMTBIT]  =  ONECHXMT(F) 

%  Creates  a  random  transmit  (XMTBIT)  bit  of  one  or  zero. 
%  Based  on  transmit  bit  it  returns  a  signal 
%  (S)  one  period  in  length  that  is  inphase 

%  if  XMTBIT  is  a  zero  and  1 80  degrees  out  of 
%  phase  if  xmtbit  is  one. 

%  decide  which  bit  transmitter  will  simulate 
%  signal  to  send 

T=1/f;  %  compute  period  of  signal 

t=linspace(0,T);  %  generate  100  time  intervals 
x=  rand  -  0.5; 


if  x  <=  0 
xmtbit=0; 

S=phase*cos(2*pi*f.*t); 

%  compute  S  for  XMTBIT  =  0 

else 

xmtbit=1 ; 

S=-1  *phase*cos(2*pi*f  .*t); 

%  compute  S  for  XMTBIT  =  1 
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end 


%  Naval  Postgraduate  School,  ECE 
%  Created  by:  LT  Paul  W.  Wehner 


3.  MACH-ZENDER  COUPLER 


function  [a,b]  =  mzcoupsm(S,T,Td,SM1  ,SM2,SM3,L) 

%  MACH  ZENDER  COUPLER  with  varied  scattering  matrices  output 
%  For  a  given  input  signal,  S  (with  a  period  T, 

%  three  periods  in  length),  returns  two  signals, 

%  A  and  B  (both  one  period  in  length)  where 
%  A  =  S(t-T_actual)+S(t)  and  B=S(t-T_actual)-S(t). 

%  T_ACTUAL  indicates  the  actual  length 
%  (i.e., seconds)  of  the  coupler  delay 

%  element  and  SM#  indicates  the 
%  scattering  matrix  of  the  individual  #'rd 
%  2x2  optical  couplers.  L  is  the  attenuation 
%  in  channel  a  to  equalizer  signal  strength 
%  lost  in  channel  b  due  to  the  last  phase 
%  shift. 

%  determine  #  discrete  quantities  in 
%  a  signal  period  (T)  of  S 
n  =  length  (S)/3; 

%  demux  signal  S  in  1st  2x2  coupler 
[al  ,b1]  =  optcoup2(SM1  ,S); 

%  delay  signal  in  channel  a 
[a2,d]  =  delay(a1  ,T,Td); 

%  adjust  length  of  channel  a  and 
%  channel  b  signals 
a3  =  a2(d-n+1  :d); 
b2  =  b1(d-n+1:d); 

%  2nd  2x2  coupler  output 
[a4,b3]  =  optcoup2(SM2,a3,b2); 

%  phase  shft  element  output 
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[q,b4]  =  optcoup2(SM3,b3); 

%  attenuate  channel  a 
[a5]  =  atten(a4,L); 

%  define  MZ  coupler  output 
a  =  real(a5); 
b  =  real(b4); 

%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


4.  OPTICAL  COUPLER 


function  [a,b]  =  optcoup2(SM,S1,S2) 

%  optical  2x2  coupler 
%  OPTCOUP2(SM,S1  ,S2) 

%  For  a  given  input  signal,  SI ,  and  a  possible 
%  second  signal  S2,  provides  two  outputs, 

%  A  and  B.  Both  outputs  are  attenuated  by  a 
%  factor  of  3dB  and  multiplied  by  the  scattering 
%  matrix,  SM.  SM  must  be  [2x2]. 

%  determime  if  second  signal  S2  was  input 
if  nargin  ==  2 
S2=zeros(1  ,length(S1 )); 
elseif  length  (SI)  ~=  length  (S2) 
disp('ERROR  -  Input  arrays  (SI  and  S2)  must') 
disp(’  be  of  the  same  length') 
a=0; 
b=0; 
return; 
end 

%  define  input  matrix 
y  =  [SI  ;S2]; 

%  multiply  by  scattering  matrix 
x  =  SM*y; 


%  attenuate  output  by  3dB 
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x  =  x./sqrt(2); 


%  define  individual  outputs 
a  =  x(1 
b  =  x(2,:); 


%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


5.  DELAY  ELEMENT 


function  [a,d]  =  delay  (a,T,Td) 

%  DELAY  ELEMENT 

%  DELAY(A,T,Td),  for  a  given  signal  A  (A  must 
%  be  three  periods  in  length  with  a  period  of 
%  T  seconds)  outputs  a  new  signal  A  delayed  Td 
%  seconds.  Indicates  at  which  element 
%  value,  N,  in  a  the  delayed  vector  begins 
%  where  the  vector  is  listed  right  to  left 

%  create  working  copy  of  a 
b  =  a; 

%  determine  #  discrete  quantities  of  one 
%  period  of  a 
i  =  length  (a)/3; 

%  compute  delay  as  percentage  of  actual  delay 
perc_Td  =  Td/T; 

%  equate  percent  delay  to  #  of  discrete 
%  quantities  of  a 

j  =  round(i  *  perc_Td);  %  quantized  delay 
k  =  (3*i)  -  j;  %  new  start  of  a 


%  zero  out  delayed  discrete  elements  of  a 
for  f  =  (3*i):-1:(k+1) 
a(l)=0; 
end 
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m  =  3*i; 


%  overwrite  a  with  delayed  version  of  itself 
for  n=  k:-1  :(2*i-j+1 ) 
a(n)=b(m); 
m=m-1 ; 
end 

d=k; 

%  Naval  Postgraduate  School,  ECE 
%  Created  by:  LT  Paul  W.  Wehner 


6.  ATTENUATOR 

function  [a]  =  atten(a,L) 

%  ATTENUATOR 
%  ATTEN(A,L) 

%  For  given  input  signal,  A,  performs 
%  attenuation  on  signal  equal  to 
%  L  dB. 

a  =  a  /  (sqrt(1 O'XL/I 0))); 

%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


7.  ELECTRICAL  CIRCUIT  ELEMENTS 

function  [rcvbit]  =  rcvr(xa,xb,R,Zo,No,T) 

%  OPTICAL  DPSK  SIGNAL  RECEIVER  OUTPUT 
%  [RCVBIT]=CRCVR(XA,XB,R,ZO,NO,T) 

%  For  given  input  signals,  XA  and  XB 
%  this  applies  photodetection 
%  proportional  to  the  responsivity  (R)  to 
%  convert  the  optical  signal  to  electrical, 

%  amplifies  the  electrical  signal 
%  with  a  gain  of  (Zo)  and  then 
%  integrates  the  signal  over  a 
%  single  time  period  (T)  to  obtain  signal 
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%  energy  in  each  branch.  A  comparison  of 
%  two  branches  is  then  made  to 
%  determine  the  transmitted  bit.  Noise  is 
%  added  in  the  photodetector  and 
%  the  amplifier  with  a  power  spectral 
%  density  of  No/2. 

%  inut  optical  signal  to  photodetector 
xal  =  pd(xa,R,No,T); 
xbl  =  pd(xb,R,No,T); 

%  input  signal  to  transimpedance  amplifier 
xa2  =  tranamp(xa1  ,Zo,No,T); 
xb2  =  tranamp(xb1  ,Zo,No,T); 

%  integrate  both  channels  over  T 
xa3  =  intg(xa2,T); 
xb3  =  intg(xb2,T); 

%  compare  channels  to  decide  received  bit 
if  xa3  >  xb3 
rcvbit  =  1 ; 
else 

rcvbit  =  0; 
end 

%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


8.  PHOTODETECTOR 

function  [a]  =  pd(x,R,No,T) 

%  PHOTODETECTOR 
%  PD(x,R,No,T) 

%  For  given  input  signal,  X,  performs 
%  photdetection  on  signal  power  proportional 
%  to  the  responsivity,  R,  in  amps/watt. 

%  Assume  signal  is  output  over  a  1  ohm 
%  load  such  that  amps  correlate  directly  to 
%  volts.  Adds  photodetection  shot  noise  to 
%  output  signal  A  with  PSD(No/2)  as  defined 
%  from  user  input. 
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%  Create  shot  noise 
s  =  randn(1  ,length(x)); 

%  adjust  shot  noise  PSD 
amp  =  sqrt(No/T); 
s  =  amp  *  s; 

%  compute  signal  power 
Px  =  (x.A2)./2; 


%  multiply  input  power  by  responsivity 
%  to  get  output  in  amps 
Ax  =  R  *  Px; 

%  add  shot  noise  to  create  output  signal 
a  =  Ax  +  s; 

%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


9.  TRANSIMPEDANCE  AMPLIFIER 

function  [a]  =  tranamp(x,Zo,No,T) 

%  TRANSIMPEDANCE  AMPLIFIER 
%  TRANAMP(x,Zo,No,T) 

%  For  given  input  signal,  X,  performs 
%  uniform  amplification  on  signal 
%  with  gain  Zo.  Adds  noise  to  signal  prior 
%  to  amplification  with  PSD  No/2. 

%  Assume  signal  is  output  over  a  1  ohm 
%  load  such  that  amps  correlate  directly  to 
%  volts.  Adds  photodetection  shot  noise  to 
%  output  signal  A  with  PSD(No)  as  defined 
%  from  user  input. 

%  Create  amp  noise 
s  =  randn(1  ,length(x)); 


%  adjust  amp  noise 
amp  =  sqrt(No/T); 
s  =  amp  *  s; 
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%  add  amp  noise  to  input  signal 
a  =  x  +  s; 

%  adjust  output  by  gain,  Zo 
a  =  Zo  *  a; 


%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


10.  INTEGRATOR 

function  [a]  =  intg(x,T) 

%  INTEGRATOR 
%  INTG(x,T) 

%  For  given  input  signal,  X,  performs 
%  trapazoidal  integration  on  signal 
%  over  one  time  period  T.  Assumes 
%  elements  of  X  are  equally  spaced 
%  from  0  to  T. 

a  =  (T/length(x))  *trapz(x); 


%  Naval  Postgraduate  School 
%  Created  by:  LT  Paul  W.  Wehner 


11.  Q  FUNCTION 

function  [Q,P]  =  qfn(x) 

%  This  function  computes  the  Gaussian  Q-function 
%  using  the  rational  approximation  26.2. 17  of 
%  Abromowitz  and  Stegun.  For  sufficiently  large 
%  arguments,  the  asymptotic  expansion  26.2. 1 2  is  used. 
%  The  rational  approximation  is  accurate  to  within 
%  7.5e-8.  Q  is  the  value  of  the  Q-function  (integral  of 

%  Gaussian  pdf)  of  zero  mean  and  unit  variance  from  the 
%  input  quantity  to  infinity,  and  P  =  1  -  Q  is  the 
%  emulative  distribution  function.  The  call  is 
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%  [Q,P]  =  qfn(x). 

% 

%  R.E.  Ziemer 
%  4/17/94 

b(l)=0.3 1938 153; 
b(2)=-0.356563782; 
b(3)=l. 78 1477937; 
b(4)=-l. 821255978; 
b(5)=l. 330274429; 
p=0.2316419; 
y=abs(x); 

T=l/(l+p*y); 

TT=zeros(l,5); 
for  r=l:5 
TT(r)=TAr; 
end 

Z=exp(-yA2/2)/sqrt(2*pi); 
if  y<4 

Ql=Z*b*TT'; 

else 

Q 1  =(Z/y)*(  1  -  l/yA2+3/yA4); 
end 
if  x>0 

Q=Q1; 

else 

Q=i-Qi; 

end 

P=l-Q; 
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